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Abstract 

The multivariate resultant is a fundamental tool of computational 
algebraic geometry. It can in particular be used to decide whether a 
system of n homogeneous equations in n variables is satisfiable (the 
resultant is a pol5momial in the system's coefficients which vanishes 
if and only if the system is satisfiable). In this paper, we investigate 
the complexity of computing the multivariate resultant. 

First, we study the complexity of testing the multivariate resultant 
for zero. Our main result is that this problem is NP-hard under de- 
terministic reductions in any characteristic, for systems of low-degree 
pol5momials with coefficients in the ground field (rather than in an 
extension). In null characteristic, we observe that this problem is in 
the Arthur-Merlin class AM if the generalized Riemann hypothesis 
holds true, while the best known upper bound in positive character- 
istic remains PS PACE. 

Second, we study the classical algorithms to compute the resul- 
tant. They usually rely on the computation of the determinant of an 
exponential-size matrix, known as Macaulay matrix. We show that 
this matrix belongs to a class of succinctly representable matrices, for 
which testing the determinant for zero is proved PSPACE-complete. 
This means that improving Canny's PS PACE upper bound requires 
either to look at the fine structure of the Macaulay matrix to find an 
ad hoc algorithm for computing its determinant, or to use altogether 
different techniques. 



*This work was partially funded by the European Community (7th PCRD Contract: 
PIOF-GA-2009-236197). 
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I Introduction 



Given two univariate polynomials, their Sylvester matrix is a matrix built 
on the coefficients of the polynomials which is singular iff the polynomials 
have a common root. The determinant of the Sylvester matrix is known 
as the resultant of the polynomials. This determinant is easy to compute 
since the size of the Sylvester matrix is the sum of the degrees of the poly- 
nomials. The study of the possible generalizations to multivariate systems 
comes within the scope of elimination theory [49} [38} [15} [47} [i6l. This 
theory proves that the only case where a unique polynomial can testify 
to the existence of a common root to the system is the case of n homo- 
geneous polynomials in n variables: the resultant of a square system of 
homogeneous polynomials fi, . . . ,fn G . . . ,x„] is a polynomial in 

the indeterminate coefficients of fi, . . . , f„ which vanishes if and only if 
f\, ■ ■ ■ , fn have a nonzero common root in the algebraic closure of K. The 
resultant of such a system is known as the multivariate resultant in the lit- 
erature. This captures the case of two univariate polynomials via their 
homogenization. Furthermore, in many cases a system of more than n 
homogeneous polynomials in n variables can be reduced to a system of 
n homogeneous polynomials, and so the square case is an important one. 
This result is sometimes known as Bertini's theorem. (As explained toward 
the end of this section, we will use an effective version of this result in one 
of our NP-hardness proofs.) In this paper, we focus on the multivariate 
resultant which we simply refer to as the resultant. 

The resultant has been extensively used to solve polynomial systems 
1 36 [43lfiolli2l1 and for the elimination of quantifiers in algebraically or real- 
closed fields 144} [24I . More recently, the multivariate resultant has been of 
interest in pure and applied domains. For instance, the problem of robot 
motion planning is closely related to the multivariate resultant HI [TT| , 
and more generally the multivariate resultant is used in real algebraic 
geometry [9, 29]. Finally, in the domain of symbolic computation progress 
has been made for finding explicit formulations for the resultant [28, [14} 
[3[i2l[25|, see also 

The aim of this paper is to contribute to the study of the complexity 
of the resultant. We both study the complexity of testing the resultant for 
zero, that is the satisfiability of a system of n homogeneous polynomials 
in n variables, and the complexity of explicitly (and exactly) computing 
the resultant. 

A shorter version of this paper [22] has been published in the Proceed- 
ings of MFCS 2010. It contains material from Sections [2] and [3] Material 
of Section [4] appeared in preliminary form (and in French) in |2T|, but was 
not published so far. 
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1.1 Definitions and notations 



In this paper, K denotes any field and ]fC[a:i, . . . , the ring of polynomials 
in n indeterminates over K. The algebraic closure of K is denoted by K. 
For a prime number p, the finite field with p elements is denoted by Fp. 
The notation is extended to the characteristic zero by Fq = Q. 

A tuple (fli,...,fls) is denoted by a when there is no possible con- 
fusion on the range of the index. In particular, for a tuple of integers 
a = («!,...,«„), x'^ denotes the monomial x'^^x'^ ■ ■ ■ xfi", and in the next 
definition, 7 is the tuple of the 7;,a's. The total degree of the monomial x"^ 
is |a| = ai + • ■ • + DCn- 

Definition 1. Let K be a field and /i, . . . ,/n be n homogeneous polyno- 
mials in ]K[:ci, . . . , fi{x) = '}2\oc\=di liA^"^- ^he multivariate resultant R of 
/i, . . . ,/„ is an irreducible polynomial in K[j] such that 

R(7) =0 ^ 3xe K'\x ^ OJi{x) = ■■■= fn{x) = 0. (1) 

The multivariate resultant is unique up to a constant factor. 

The existence of the multivariate resultant is not evident. We refer 
to [49} [35! for a proof of this fact. The uniqueness comes from the fact 
that two irreducible polynomials having the same roots are equal up to a 
constant factor. 

The first problem we are interested in is testing the resultant for zero. 
This is the same as deciding whether a square system of homogeneous 
polynomials (that is n polynomials in n variables) has a non-trivial root. 
This is closely related to the decision problem for the existential theory of 
an algebraically closed field. This problem is sometimes called the Hilbert 
Nullstellensatz problem: 

Definition 2. Let K be a field. The Hilbert Nullstellensatz problem over K, 
HN(K), is the following: Given a system / of s polynomials in K[xo, . . . , 

— n+l 

does there exist a root of / in K ? 

Let us now assume that the s components of / are homogeneous poly- 
nomials. Then the homogeneous Hilbert Nullstellensatz problem over K., 
H2N(]K), is to decide whether a non trivial (that is, nonzero) root exists in 

If / is supposed to contain as many homogeneous polynomials as vari- 
ables, the problem is called the Resultant over K, RESULTANT(]fC). 

In the case of the field Q, it is more natural to have coefficients in Z. 
We shall use the notations HN, H2N and Resultant for the case where 
the system is made of integer polynomials. 

In the case of polynomials with coefficients in Z, Canny [8] gave in 1987 
a PS PACE algorithm to compute the resultant. To the authors' knowledge. 
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this is the best known upper bound. In this paper we show that testing 
the resultant for zero is NP-hard in any characteristic. In other words, 
Resultant(]K) is NP-hard for any field K. 

Main results and proof techniques 

Sections E] and [3] are devoted to the study of the decision problem, first in 
null characteristic and then in any characteristic. Section [4] focuses on the 
evaluation problem. 

In Section [27T1 we observe that for polynomials with integer coefficients, 
testing the resultant for zero is a problem in the Arthur-Merlin (AM) class. 
This result assumes the Generalized Riemann Hypothesis (GRH), and fol- 
lows from a simple reduction to the Hilbert Nullstellensatz. For this prob- 
lem, membership in AM assuming GRH was established in [30]. 

In characteristic zero, it seems to be a "folklore" result that testing the 
resultant for zero is NP-hard, and one can find a proof of this fact in p3| . 
A similar result can be obtained by considering a system of two homo- 
geneous polynomials, but given in lacunary representation (their degree 
can therefore be exponential in the input size). This result of incompara- 
ble strength is a reformulation of a theorem of Plaisted p^ . We give the 
proofs of these two results to be able to argue about their irrelevance in 
positive characteristic. 

The first proof does not carry over to positive characteristic since it is 
a reduction from the problem Partition [[20I Problem SP12] whose NP- 
hardness relies in an essential way on the fact that the data are integers (in 
fact, in any finite field the analogue problem can be solved in polynomial 
time by dynamic programming). 

Plaisted's result can be adapted to positive characteristic |j5o} EH but 
this requires randomization. By contrast, our ultimate goal is N P-hardness 
for deterministic reductions and low degree polynomials. We therefore 
need to use different techniques. Our starting point is a fairly standard 
encoding of 3-SAT by systems of polynomial equations. Using this encod- 
ing we show at the beginning of Section [3] that deciding the existence of 
a nontrivial solution to a system of homogeneous equations is NP-hard 
in any characteristic. The resulting system has in general more equations 
than variables. In order to obtain a square system we explore these two 
basic strategies: 

(i) Decrease the number of equations. 

(ii) Increase the number of variables. 

In Section [3TT] we give a randomized N P-hardness result based on the first 
strategy. The idea is to replace the initial system by a random linear combi- 
nations of the system's equations (the fact this does not change the solution 
set is sometimes called a "Bertini theorem"). 



4 



In Section [3T2I we use the second strategy to obtain two N P-hardness 
results for deterministic reductions. The main difficulty is to make sure 
that the introduction of new variables does not create spurious solutions 
(we do not want to turn an unsatisfiable system into a satisfiable system). 
Our solution to this problem can be viewed as a derandomization result. 
Indeed, it can be shown that the coefficients of the monomials where the 
new variables occur could be chosen at random. It would be interesting to 
find out whether the proof based on the first strategy can also be deran- 
domized. 

Section [4] is devoted to the study of the Macaulay matrices which are 
used in several formulations of the resultant. These matrices have an 
exponential-size (in the number of variables and the degrees of the poly- 
nomials) but can be efficiently represented by a polynomial-size boolean 
circuit. This notion of circuit representation goes back to fTg\ and was de- 
veloped for graphs. We prove that computing the determinant (or testing 
it for zero) of a matrix which is given by a circuit is PSPACE-complete 
in any characteristic. This shows that to improve the best known upper 
bound on the computation of the resultant in positive characteristic, one 
needs either to use the fine structure of the Macaulay matrices, or to find 
an altogether different method. So far, this has been achieved only in 
characteristic [30] . 

2 Complexity of the resultant in null characteristic 
2.1 The resultant lies in AM 

In this section we show that testing the resultant for zero is reducible to 
HN(]K). In the case K = Z, this allows us to conclude (under the Gener- 
alized Riemann Hypothesis) that our problem is in the polynomial hierar- 
chy, and more precisely in the Arthur-Merlin class. In fact, we show that 
this applies more generally to the satisfiability problem for homogeneous 
systems (recall that testing the resultant for zero corresponds to the square 
case). 

Proposition 1. For any field K, the problem ¥iJ^{K) is polynomial-time many- 
one reducible to HN(IK). 

Proof. Consider an instance / of H2N(K), that is s homogeneous polyno- 
mials /i, . . . ,/s G K[xi, . . . ,Xn]. The polynomials fi,...,fs can be viewed 
as elements of K[xi, . . . ,x„,yi, . . . , y„] where 1/1, . . . , y„ are new variables 
which do not appear in the f/s. Let g be the system containing all the 
fi's and the new (non-homogeneous) polynomial Ylf^i Xiyi — 1. This is an 
instance of the problem HN(]K). It remains to prove that / and g are 
equivalent. 
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Given a root {ai, . . . ,a„,bi, . . . ,bn) of g, the new polynomial ensures 
that there is at least one nonzero a,. So (fli, . . . , fl„) is a non trivial root of 
/. Conversely, suppose that / has a non trivial root (fli, . . . , fl„), and let i 
be such that a, 7^ 0. Then the tuple (fli, . . . ,a„,0, . . . ,0, a^^, 0, . . . , 0) where 
aj^ corresponds to the variable j/, is a root of g. 

Thus H2N(K.) is polynomial-time many-one reducible to HN(]K). □ 

Koiran [30] proved that HN G AM under the Generalized Riemann 
Hypothesis. We denote here by AM the Arthur-Merlin class, defined by 
interactive proofs with public coins (see [2]). Thereby, 

Corollary 1. Under the Generalized Riemann Hypothesis, H2N is in the class 
AM. 

In positive characteristic, the best upper bound on the complexity of 
the Hilbert Nullstellensatz known to this day remains PS PACE (in partic- 
ular it is not known whether the problem lies in the polynomial hierarchy, 
even assuming some plausible number-theoretic conjecture such as the 
generalized Riemann hypothesis). 

2.2 The resultant is NP-hard in characteristic 

We now give a first N P-hardness result, for the satisfiability of square sys- 
tems of homogeneous polynomial equations. The first part of the theorem 
can be easily deduced from [23, Proposition 10]. The second part shows 
that the problems remains NP-hard even for systems with small integer 
coefficients (i.e., coefficients bounded by 2). This is achieved by a standard 
trick: we introduce new variables in order to "simulate" large integers 
coefficients. It is interesting to note, however, that a similar trick for re- 
ducing degrees does not seem to apply to the resultant problem (more on 
this after Theorem [2]|. 

Theorem 1. The problem Resultant of deciding whether a square system of 
homogeneous polynomials with coefficients in Z has a non trivial root is NP-hard. 

The problem remains HP -hard even if no polynomial has degree greater that 2 
and even if the coefficients are hounded by 2. 

Proof. The reduction is done from Partition which is known to be NP- 
hard [[201 Problem SP12]: Given a finite set S and a non negative integer 
weight w{s) for each s G S, the problem is to decide the existence of subset 
S' such that YLs^^s' ^i'(s) = Yosts' '^{s)- That is, the aim is to cut A into two 
subsets of same weights. 

Given such an instance of Partition where S = {si,...,s„}, let us 
define a system of polynomials. For \ <i < n, fi{x) = Xq — xf. And 

fo{xo,...,Xn) = w{Si)Xi+w{S2)x2 H \-w{Sn)x„. (2) 
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A tuple (flo/ ■ ■ ■ / fln) has to verify a,- = ±flo for each i to be a solution, hence 
the only case to consider is aq = 1 arid fl, = ±1 for i > 1. Then it is clear 
that the system has a non trivial solution if and only if S may be split into 
two subsets of equal weights. 

For the second part of the theorem, it remains to show that the coeffi- 
cients in the system can be bounded by 2. As the a;(s,)'s may be large inte- 
gers, they have to be replaced by variables. Let us write w(s,) = EyLg WjjlL 
For each iVij, a new variable W,y is introduced. For every i, the values of 
the Wij's are defined by a descending recurrence: 



Then /o is replaced by Wi^o^i + ^2,0^2 + • ■ • + Wn,oXn- Doing so, the num- 
ber of polynomials remains the same as the number of variables. Hence, 
this algorithm build a new homogeneous system where the polynomials 
have their coefficients bounded by 2 and their degrees too. One can read- 
ily check that the new system has a non trivial solution if and only if the 
original one has. In particular, if Xq is set to zero, then all other variables 
have to be set to zero too. □ 

A related result is Plaisted's [42] on the N P-hardness of deciding whether 
the gcd of two sparse univariate polynomials has degree greater than one. 
By homogenization of the polynomials, this is the same problem as in 
Theorem [T] for only two bivariate polynomials. Note that the polynomials 
are sparse and can be of very high degree since exponents are written in 
binary (this polynomial representation is sometimes called "supersparse" 
or "lacunary" E^ll. If both polynomials were dense, the resultant could 
be computed in polynomial time since it is equal to the determinant of 
their Sylvester matrix. Plaisted's theorem stated in the same language as 
Theorem [1] is the following: 

Theorem 2. Given two sparse homogeneous polynomials in Z[x, y\, it is HP -hard 
to decide whether they share a common root in 

We briefly sketch Plaisted's reduction since it will help understand 
the discussion at the end of this section. For a full proof (including a 
correctness proof), see [ I42} Theorem 5.1]. 

Proof sketch. The idea is to turn a 3-SAT instance into a system of two 
univariate polynomials which share a common root if and only if the 3- 
CNF formula is satisfiable. 

To every variable Xj is associated a prime pj, and let M = Y[j Pj where 
the product ranges over all the variables that appear in the formula. A 




(3) 
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formula (p is turned into a polynomial according to the following rules. 
A non negated variable Xy is turned into Pxji^) = x^^fi — 1 and a negated 

variable -^X^ into P-.x,(^) = 1 + x^^P' H h x^P^-^^^^P':. Then a formula 

(pV ip is turned into P<^vi/> = lcm(P^, P^). A conjunction (p = /\j (pi is turned 
into the polynomial 

P^{x) = x^J2P'P.i^W/^) (4) 

i 

This defines the first polynomial P. The second polynomial is simply 
y.M _ -[^ jYxe proof that those two polynomials share a common root if and 
only if is satisfiable is omitted. 

To obtain the result in the way we stated it, it is sufficient to homoge- 
nize P{x) and x^ — 1 with the second variable y. □ 

Theorems [T] and E] seem to be incomparable. In particular, it is not 
clear how to derive Theorem [i] from Theorem [H A natural idea would be 
to introduce new variables and use the repeated squaring trick to reduce 
the degrees of the polynomials occurring in Plaisted's result. However, as 
we now explain this can lead to the creation of unwanted roots at infinity 
in the resulting polynomial system. 

Assume for instance that we wish to get rid of all occurrences of x^ in 
a polynomial. One can add a new variable X2, replace the occurrences of 

by X2 and add a new polynomial X2 — x^. In order to keep the system 
homogeneous, the idea is to homogenize the latter polynomial: xoX2 — x^. 
The problem with this technique is that it adds some new roots with all 
variables but X2 set to 0, and in particular the homogenization variable xq. 

To give an explicit example of the problem mentioned above, let us 
consider the formula 

(xvy)A(-x)A(-y). (5) 

Let us associate the prime number 2 to the variable X, and 3 to Y (M in 
the previous proof is therefore 6). By Plaisted's construction, X is turned 
into _ i — ^■3 _ y ^^^.^ ^2 _ ^ xheir negations -iX and -lY are 

respectively turned into 1 + x^ and \ -\- x^ -\- x'^ . The disjunction of X and 
y is turned into the 1cm of x'' — 1 and x^ — 1, that is (x^ — l)(x^ + x + 1). 
Finally, we have to apply formula (j4| with the latter polynomial, 1 + x'^ and 
\ + x^ + x'^. Therefore, the two polynomials of Plaisted's construction are 
y-M _i _y3 j^x'^^ 2x^ + 9x^ + + x^- x^. It can be checked 

that as expected, those two polynomials do not share any common root. 

Applying the repeated squaring trick with homogenization on this ex- 
ample gives the following system where the two first pol}momials repre- 
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sent the original ones and the other ones are new ones: 

-X3 + X4 + 2X5 + 9xs + 2x7 + xs - = 
xe- xo = 0; X0X2 - x^ = 0; xqx^ - X2X = 
Xo^4 — = 0; X0X5 — x^x = 0; xoX(, — = 
^ X0X7 - X4X3 = 0; xqXs - xl = 0; xqx^ - xgx = 

But in that example, one can easily check that solutions with xq = ex- 
ist. Namely if we set xg and xg to the same nonzero value and all other 
variables to 0, this defines a solution to the system. 

To the authors' knowledge, there is no solution to avoid these un- 
wanted roots. Furthermore, Plaisted's result works well with fields of 
characteristic 0, but as it uses the fact that a sum of non negative terms is 
zero if and only if every term is zero, this generalizes not so well to pos- 
itive characteristic. In particular, generalizations to positive characteristic 
require randomization (see [26 J and [50]). By contrast, two of the reduc- 
tions given in the next section are deterministic and they yield systems 
with polynomials of low degree (i.e., of linear or even constant degree). 



3 NP-hardness in arbitrary characteristic 

In this section we give three increasingly stronger N P-hardness results for 
testing the resultant. As explained in the introduction, we first provide in 
Section [3TT] a NP-hardness proof for randomized reductions. We then give 
in Section [3T2I two N P-hardness results for deterministic reductions: the 
first one applies to systems with coefficients in an extension of the ground 
field, and the second (stronger) result to systems with coefficients in the 
ground field only. The starting point for these three NP-hardness results 
is the following easy lemma. 

Lemma 1. Given a field K of any characteristic, it is N P-hard to decide whether 
a system of s homogeneous polynomials in ]K[xo, . .. ,Xn] has a non trivial root. 
That is, H2N(]K) is NP-hard. 

In [ [321 , this lemma is proved using a reduction from the language 
BooLSYS. An input of Boolsys is a system of boolean equations in the vari- 
ables Xi, . . . ,X„ where each equation is of the form X, = True, X, = -iXy, 
or Xi = Xj V X/f. The question is the existence of a valid assignment for the 
system, that is an assignment of the variables such that each equation is 
satisfied. This problem is easily shown N P-hard by reduction from 3-SAT. 
We now give a proof of this lemma since the specific form of the systems 
that we construct in the reduction will be useful in the sequel. This proof 
is a slight variation on the proof from I32J. 

Proof. Let K be a field of any characteristic p, p being either zero or a 
prime number. At first, p is supposed to be different from 2. The proof 
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has to be slightly changed in the case p = 2 and this case is explained at 
the end of the proof. 

Let B be an instance of Boolsys. Let us define a system of homoge- 
neous polynomials from this instance with the property that B is satisfiable 
if and only if the polynomial system has a non trivial common root. The 
variables in the system are Xq, . . . ,x„ where x,-, 1 < i < n, corresponds to 
the boolean variable X, in Boolsys, and Xq is a new variable. The system 
contains four kinds of polynomials: 

• Xq — xf, for each i > 0; 

• Xq ■ {xi + Xq), for each equation = True; 

• Xq ■ {xi + Xj), for each equation Xi = -iXy; 

• (xi + xo)'^ — {xj + Xq) ■ + Xq), for each equation X, = Xy V X]^. 

Let us denote by / the polynomial system obtained from B. The first 
kind of polynomials ensures that if {uq, . . . ,a„) is a non trivial root of /, 
then Aq = flf = • • • = a^. Now if / has a non trivial root {uq,. . . ,an), 
then one can readily check that the assignment X/ = True if fl; = —uq and 
X, = false if a, = uq satisfies B. Conversely, if there is a valid assigrunent 
Xi, . . . , X„ for B, any (n + l)-tuple {uq, . . . ,an) where flo 7^ and ai = —Uq 
if Xj = True and fl, = uq if X, = false is a non trivial root of /. 

This proof works for any field of characteristic different from 2. The 
problem in characteristic 2 is the implementation of Boolsys in terms of 
a system of polynomials. Indeed, for the other characteristics, the truth 
is represented by —uq and the falseness by uq. In characteristic 2, those 
values are equal. Yet, one can just change the polynomials and define in 
the case of characteristic 2 the following system: 

• xoXi — xf, for each i > 0; 

• xo{xi + Xq), for each equation X, = True; 

• Xo(x; + Xj + Xq), for each equation X, = -iXy; 

• xf + XjXk + Xq • {xj + Xk), for each equation X, = Xj V X^. 

Now, given any nonzero value aq for xq, the truth of a variable X,- is rep- 
resented by Xi = flo whence the falseness is represented by x, = 0. A root 
of the system is in particular a root of the polynomials defined by the first 
item. Therefore each has to be set either to aq or to 0. The system has a 
non trivial root if and only if the instance of Boolsys is satisfiable. □ 
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3.1 A randomized reduction 

We now give the first of our three N P-completeness results in positive 
characteristic. The proof also applies to characteristic zero, but in this 
case Theorem [T] is preferable (its proof is simpler and the N P-hardness 
result stronger since it relies on deterministic reductions). For more on 
randomized reductions, see [2]. 

We begin with a result on algebraic varieties in algebraic closed fields. 
This is a classical result in algebraic geometry see also fp^ Proposi- 
tion 1] for a proof in our language. 

Theorem 3. Let K be an algebraically closed field and V an algebraic vari- 
ety o/]K"+^ defined by a set of homogeneous degree-d polynomials fi,...,fs S 
K[x(j, . . . , x„]. This variety can be defined by (n + 1) homogeneous degree-d poly- 
nomials gi, ...,gn+i G IK[xo, . . . , Xn]. Moreover, suitable gi's can be obtained by 
taking generic linear combinations of the fi. That is, we can take gi = J^j^i O'-ijfj 
where [otij) is a matrix of elements o/K., and the set of suitable matrices is Zariski- 
dense in 

This result leads us to out first N P-hardness result in arbitrary charac- 
teristic. 

Theorem 4. Let p be either zero or a prime number. The following problem is 
N P-hard under randomized reductions: 

• Input: A square system of homogeneous equations with coefficients in a 
finite extension o/Fp. 

• Question: Is the system satisfiable in the algebraic closure of¥p? 

In the case p = 0, the results also holds for systems with coefficients in Z. 

Proof Lemma [1] shows that it is N P-hard to decide whether a non square 
polynomial system / with coefficients in Fp has a non trivial root. From /, 
a square system g is built in randomized polynomial time. 

Let us denote by fj, 1 < j < s, the components of /. They are homo- 
geneous polynomials in Fp[xo/ ■ ■ - fXn]. The components of g are defined 
by 

s 

for < / < n. In the sequel, we explain how to choose the a,y's for / and g 
to be equivalent. For any choice of the ocifs a root of / is a root of g. Thus 
it is sufficient to show how to choose them so that g has no non trivial root 
if the same is true for /. 
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The property the a^'s have to satisfy is expressed by the first-order 
formula 



n s 



^ Vxo • • ■ Vx„ 1^ A /;(^) = j V V E ^ Oj • (8) 

The formula O belongs to the language of the first-order theory of the 
algebraically closed field F^. This theory eliminates quantifiers and 
is therefore equivalent to a quantifier-free formula of the form 

Y(^^) = \J {/\P^,{^)=Q^f\ Q,„(«) ^ ) , (9) 
\ ; m / 

where P^i, Qkm G W] ■ As a special case of [iF. Theorem 2], one can bound 
the number of polynomials in Y as well as their degrees by 2P°'y("'i°8(^+")) 
where poly represents some polynomial independent from O. 

Theorem [3] shows that the set A of tuples satisfying O is Zariski-dense 

in Fp . Since A is dense, and A is also defined by Y, one of the clauses 

of Y must define a Zariski dense subset of Fp . This clause is of the 
form Am Qm(a) ^ 0. 

To satisfy O, it is sufficient for the cty's to avoid the roots of a polyno- 
mial Q = Ylm Qm- As mentioned before, it is known that Y contains at 
most 2P°iy("'i°g(^^+")) polynomials of degree at most 2P°iy("'iog(s+«)). Thus, Q 
is a polynomial of degree at most 2^P°'y("''°8(^+")). Consider now a finite 
extension K of Fp with at least 22+2p°iy(«.iog(s+n)) elements (that is, of poly- 
nomial degree). If we choose the a,y's uniformly at random in K, then with 
probability at least 3/4 they are not a root of Q (by the Schwartz-Zippel 
Lemma). Thus with the same probability, they satisfy O. Note that IK can 
be built in polynomial-time with Shoup's algorithm | |46| when p is prime 
(for p = 0, we take of course K = Q). 

To sum up, we build from / a square system g defined by random 
linear combinations of the components of /. If / has a non trivial root, 
then it is a root of g too. Conversely, if / has no non trivial root, then with 
probability at least 3/4 it is also the case that g has no nontrivial root. □ 

In characteristic zero the bounds in the above proof can be sharpened: 
instead of appealing to the general-purpose quantifier elimination result 
of |fi8| we can use a result of [34] . Indeed, it follows from Section 4.1 
of I34J that there exists a polynomial F of degree at most 3"+^ such that 
F(a) 7^ implies that g has no non trivial root as soon as it is true for /. 
This polynomial plays the same role as Q in the previous proof but the 
bound on its degree is sharper. 



12 



3-2 Deterministic Reductions 

We now improve the N P-hardness result of Section [3TTI we show that the 
same problem is NP-hard /or deterministic reductions. This result is not 
only stronger, but also the proof is more elementary (there is no appeal to 
effective quantifier elimination). 

Recall from the introduction that for a field K, Resultant (K) is the 
following problem: 

• Input: A square system of homogeneous equations with coefficients 
inK. 

• Question: Is the system satisfiable in the algebraic closure of K.? 

Theorem 5. Let p be either zero or a prime number. There exists a finite extension 
IK o/Fp such that Resultant(K) is HP-hard. 

In the case p = 0, this result also holds for systems with coefficients 
in Z. 

Proof. The proof of Lemma [1] gives a method to implement an instance of 
BooLSYS with a system / of s homogeneous polynomials in n + 1 variables 
with coefficients in Fp. It remains to explain how to construct a square 
system g that has a non trivial root if an only if / does. Let us denote by 
/i, . . . ,/s the components of /, with for each i = 1,. . . ,n, fi = Xq — xf if 
p 7^ 2 and /, = xqx, — xf if p = 2. A new system g of s polynomials in s 
variables is built. The s variables are xq, . . . , x„ and yi, . . . ,ys-n-i, that is 
(s — n — 1) new variables are added. The system g is the following: 

fn{x) 

f„+i{x) +Xyl 

fn+i{x) -y\ +\yl 

. (10) 

fn+iix) -y\_^ +Ay2 

M^) -Vl-n-l I 

The parameter A is to be defined later. Clearly, if / has a non trivial root 
a, then (a, 0, ... , 0) is a non trivial root of g. Let us now prove that the 
converse also holds true for some A: \i g has a non trivial root, then so 
does /. Note that a suitable A has to be found in polynomial time. 

Let (flo/ • • • / ^1/ • • • / fes-n-i) be any non trivial root of g. In particular, 
/i(fl) = ■ • ■ = fn{a) = 0. Hence Aq = • ■ • = fl^ if p 7^ 2, and fl, G {0,flo} 



g{x,y) = 
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for every / if p = 2. Now, either uq = and fi{a) =0 for every /, or aq can 
be supposed to equal 1. Therefore, if p / 2 either a = or a,- = ±1 for 
every /, and if p = 2 either fl = or fl, G {0, 1} for every i. Let us define 
£i = fn+i{^) G Fp. As {a, b) is a root of g, the &?'s satisfy the linear system 

ei + AYi = 0, 

€2 - Yi + AY2 = 0, 

; (11) 

£s-n-l ~ Ys-n-2 + -^'^s-n-l = 0, 
£s— « Xi— n— 1 ~ 0- 

This system can be homogenized by replacing each e,- by e/Yo where Yq 
is a fresh variable. This gives a square homogeneous linear system. The 
determinant of the matrix of this system is equal to 

(-1)'-"-^ (ei + e2A + • • • + es-nA^-"-') • 

Let us consider this determinant as a polynomial in A. This polynomial 
vanishes identically if and only if all the e/s are zero. In that case, the only 
solutions satisfy Y,- = for i > 0, that is (fl,0) is a root of g and therefore 
fl is a root of /. If some e,'s are nonzero, this is a nonzero polynomial of 
degree (s — n — 1). If A can be chosen such that it is not a root of this 
polynomial (for any possible nonzero value of e), then the only solution 
to the linear system is the trivial one. This mearis that the only non trivial 
root of ^ is (fl, 0) where a is a root of /. 

If the polynomials have coefficients in Z, A = 3 (or any other integer 
A > 2) satisfies the condition. Indeed, one can check that £{ = fn+i{a) G 
{—4,0,2,4} when Aq = 1- The determinant is zero if and only if e[ + 

4AH ^£'^_„A'-"-^ = where e'- = e;/2 G {-2,0,1,2}. For each i, 

let = max{e|,0} and ej^ = max{— e-,0}. Then e'- = — , and < 
< 2. Now the determinant is zero if and only if I]/e^3' = YLi^i'^^- 
By the unicity of base-3 representation, this means that for all i, ej' = , 
and so e- = 0. 

For a field of positive characteristic, this argument cannot be applied. 
The idea is to find a A that is not a root of any polynomial of degree 
(s — n — 1). Nothing else can be supposed on the polynomial because if 
p = 3 for example, any polynomial of F3 [A] can appear. This also shows 
that A cannot be found in the ground field. Suppose an extension of degree 
(s — n) is given as Fp[X]/ (P) where P is an irreducible degree-(s — n) 
polynomial with coefficients in Fp. Then a root of P in Fp[X]/(P) cannot 
be a root of a degree- (s — n — 1) polynomial with coefficients in Fp. Thus, 
if one can find such a P, taking for A the indeterminate X is sufficient. For 
any fixed characteristic p, Shoup gives a deterministic polynomial-time 
algorithm |46J that given an integer N outputs a degree-N irreducible 
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polynomial P in Fp[X]. Thus, the system g is now a square system of 
polynomials in (Fp[X]/ (P)) [x,y] and this system has a non trivial root if 
and only if / has a non trivial root. And Shoup's algorithm allows us to 
build g in polynomial time from /. 

For any field ¥p, it has been shown that from an instance B of Boolsys 
a square system g of polynomials with coefficients in an extension of Fp (in 
Z for integer polynomials) can be built in deterministic polynomial time 
such that g has a non trivial root if and only if B is satisfiable. This shows 
that the problem is NP-hard. □ 

The previous result is somewhat unsatisfactory as it requires, in the 
case of positive characteristic, to work with coefficients in an extension 
field rather than in the ground field. A way to get rid of this limitation 
is now shown. Yet, a property of the previous result is lost. Instead of 
having constant-degree (even degree-2) polynomials, our next result uses 
linear-degree polynomials. It is not clear whether the same result can be 
obtained for degree-2 polynomials (for instance, as explained at the end 
of Section [2721 reducing the degree by introducing new variables can create 
unwanted solutions at infinity). 

The basic idea behind Theorem|6]is quite simple (we put the irreducible 
polynomial used to build the extension field into the system), but some 
care is required in order to obtain an equivalent homogeneous system. 

Theorem 6. For any prime p, Resultant(Fp) is UP-hard under deterministic 
reductions. 

Proof. The idea for this result is to turn coefficient A in the previous proof 
into a variable and to add the polynomial P as a component of the system. 
Of course, considering A as a variable implies that the polynomials are 
not homogeneous anymore. Thus, it remains to explain how to keep the 
system homogeneous. 

First, the polynomial P needs to be homogenized. This is done through 
the variable xq in the canonical way. As P(A) is irreducible, it is in par- 
ticular not divisible by A. Hence, the homogenized polynomial P(A, xq) 
contains a monomial aX'^ and another one /3xq where d is the degree of P. 
Hence Xq is zero if and only if A is. 

The other polynomials have the form fn+i{x) — yj_i + Ay?. It is impos- 
sible to homogenize those polynomials by multiplying and yj_^ by xq 
(or any other variable) because then the variable y,_i never appears alone 
in a monomial, and a s-tuple with all variables set to but would be a 
non trivial solution. Moreover, in the previous proof, the fact that the y/s 
all appear with degree 2 is used to consider the system as a linear system 
in the yj. Thus replacing the monomial Ay? by Ay, does not work either. 
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Instead, we construct the slightly more complicated homogeneous system: 



fnix) 



„s— n— 1 



„s-n-2 



fn+l{x) 
'fn+2{x) 



'fn+i{x) 



„s-Ji-J+2 



-Ayr" 
-Ayr"" 



-Ay? 



,s— n— J+l 



V P(A,xo) 



'Vs-n-l 



-n-2 



-Ay; 



(12) 



Contrary to the previous proof, the y/s do not appear all at the same 
power. Yet, all the occurrences of each y, have the same degree, and we 
shall prove that this is sufficient. 

Let us prove that if / does not have any non trivial root, then neither 
does gh- Some of the observations made for g in the previous proof remain 
valid. Hence, it is sufficient to prove that a non trivial (s + l)-tuple {a, b, £) 
cannot be solution of gh whenever uq = l,b ^ and Aq = ■ • ■ = if p / 2 
or flj G {0, flo } if P = 2. By a previous remark on the polynomial P, £ can 
also be supposed to be nonzero. 

So, similarly as in the previous proof, let us define ei = al~^~'fn+i{a) e 
Fp. In the system g^, the variable y, only appears at the power {s — n — i + 
1). Therefore, given a value of a and £, the tuple {a, h,£) is a root of gii if 
and only if the fcp""'"^^'s satisfy the linear system 



^s—n—l 
^ ^s—n 



£Yi 



£Y,. 











(13) 



This is the same system as in the previous proof. Now if {£, 1) is supposed 
to be a root of P, as P is an irreducible polynomial of degree (s — n), £ 
carmot be a root of a univariate polynomial of degree less than (s — n) 
with coefficient in Fp. But the determinant of the linear system is such a 
polynomial, and thus cannot be zero. This determinant is then if and 
only if all the = 0. The same argimients as in the previous proof can be 
used to conclude that (a, b, £) can be a root of g^ if and only if a is a root 
of/. 
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Thus, from an instance B of Boolsys, a square homogeneous system 
gh of polynomials with coefficients in the ground field Fp is built in deter- 
ministic polynomial time. This system has a non trivial root if and only if 
B is satisfiable. The result is proved. □ 

4 Macaulay matrices 

There exist several formulations of the multivariate resultant. More pre- 
cisely, there are quite a lot of determinant-based formulations but they give 
in certain cases only multiple of the resultant: It can happen that the deter- 
minant that is computed vanishes even though the polynomial system has 
no root lE8l|T4} [^[T2'.'25l. Some exact formulations exist based on Macaulay 
matrices. The resultant can be computed as the gcd of the determinants 
of several Macaulay matrices (viewed as polynomials) [8J, or as the quo- 
tient of two determinants of Macaulay matrices [ j49| . Note also that this 
quotient can actually be turned into a single determinant using Strassen's 
method to eliminate divisions [27]. But in all cases, a computation of the 
determinant of one or several Macaulay matrices is needed. Canny ||8| 
gives an exact algorithm to compute the resultant based on determinants 
of Macaulay matrices that runs in polynomial space. In Valiant's algebraic 
model of computation (see [5]), the resultant polynomials belongs to the 
class VPSPACE |j33j, which can be seen as an equivalent to the boolean 
class PSPACE in this setting. 

In this section, we shall prove that an improvement of Canny's algo- 
rithm for the computation of the resultant polynomial is likely to require 
some new techniques, in particular techniques not based on Macaulay 
matrices. More precisely, we show that despite being of exponential size, 
Macaulay matrices are efficiently representable by circuits. One could thus 
wish to use this small representations to get an efficient algorithm. Yet, 
we prove that computing the determinant of such matrices is in general 
PSPACE-complete (and even testing it for zero). This holds in any char- 
acteristic. Therefore, there is an exponential blow-up of the complexity 
when the matrix is given as a circuit. Indeed, when the matrix is given in 
a standard (non succinct) way, the determinant is very closely related to 
the complexity class #L which is the counting version of L (see H for more 
on this). 

This proves that improving the PSPACE upper bound requires either 
to look at the fine structure of Macaulay matrices and use an ad hoc algo- 
rithm to compute their determinants, or to go through completely different 
methods. The complexity of computing the determinant of succinctly rep- 
resented matrices has recently been studied by Malod [40] in the settings 
of Valiant's algebraic model of computation. He shows the VPSPACE- 
completeness of this problem. 
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4.1 Representation of Macaulay matrices 

We first define the Macaulay matrices. This presentation follows [ j39| and 
||8 | . We consider a system / of n homogeneous polynomials in ]K[xi, . . . , x„] . 
Let di, d„ be the respective degrees of the /,'s, and d = 1 + Yll^i i^i — 1) 
the degree of the system. Let also Mon^ = {x* : cti -\- ■■ ■ + «.„ = d} the set of 
degree-cf monomials. Note that | Mon^f | = 

The Macaulay matrices depend on an ordering of the variables. Let 
us consider the ordering Xi ^ X2 ~< ■ ■ ■ ~< Xn- The matrix Mac is defined 
as follows: The rows and columns of Mac are indexed by the elements of 
Mon^j, ordered by the reverse lexicographic order on the tuples a: a < a' 
if there exists / such that a, < and for all j > i, ocj = oc'-. (There are 
thus I Mon^ I rows and columns.) The row of index x"' represents the 
polynomial 

— where i = min{/ : x-' divides x"^}. (14) 
X.- 

Note that {/ : x'^' divides } 7^ thanks to the definition of d. The 
other Macaulay matrices are similarly defined but with the minimum in 
the formula depending on another ordering of the variables. In particular, 
we can consider the n orderings satisfying all but one of the n inequalities 
Xi -< x;_|_i (1 < i < n) and x„ -< xi. That is, all possible ordering induced 
by the cycle (xi, . . . , x„), broken at one place. Then one can show that the 
gcd of the n Macaulay matrices obtained in this way is the resultant of the 
polynomial system [39]. 

We now aim to prove that the matrix Mac is easily representable by a 
circuit. The following definition is the straightforward adaptation of the 
notion of Small Circuit Representation for graphs |T9| , also known as succinct 
representation. Several authors have studied the transfers of complexity 
results between a (classical) graph problem and its succinct version, that is 
the same problem but where the input is a circuit describing the graph [jigj 
[4i}[37}|3|. Note that some variants of this notion have also been studied 
[4} [17), as well as a counting version | |48| and version for BSS machines HI. 

Definition 3. A circuit representation of an integer matrix M of dimensions 
(n X m) is a multiple-output circuit Cm with two inputs of [logn] and 
[logm] bits respectively, that on input (i,;) (written in binary) evaluates 
to (the binary representation of) M,y. 

A circuit representation of a graph G is a circuit representation of its 
adjacency matrix. 

Proposition 2. The matrix Mac has a circuit representation of polynomial size 
(in n and d). 

Proof. We actually give a polynomial time algorithm that on input n, d, i 
and ; outputs the entry MaCy. 



18 



The first step of the algorithm is to find the monomials x'^ and 
corresponding to the z-th row and the /-th column respectively. Let us 
write Arf = {(ai,. . . ,OLn) ■ a.i + ■ ■ ■ + a„ = d}, so that Mon^j = {x" : a G 
A^}. We need to find the z-th j-th elements of A^ in reverse lexicographic 
order. 

With the reverse lexicographic order on A,^, the first elements are tu- 
ples of the form (ai, . . . ,a„_i,0), then (ai, . . . ,a„_i, 1), and this up to 
(0, . . . , 0, d) which is the largest element. Given an index i, we first want to 
find the value of a„ of the z-th element. To this end consider, for 1 < k < n, 

= {a : o-i + ■ ■ ■ + o-k = d and a/c+i = 0}. In particu- 

lar, Arf = A^, and C for all k. Moreover, for all a G A^ and 

a' G \ A^, then a < a.'. The elements of of the form (ai, . . . , v) 
for some value v of a.„ are in bijection with A^zl- Thus, if we know that 

Efjo""^ l^r^l < ^ ^ I^to" l^r^l' we know that a„ = v. To continue, we 
remark that if the f-th element of A^ satisfies = v, then it is also the 
element of index (z — Y^i^q^^ I^"^ I ) in ^l/Z^- This allows us to find recur- 
sively all the values of the aj^'s. This gives us an algorithm, and we shall 
prove it runs in polynomial time. Since 

_ fd + k-l\ _ , d + k-1 

we can compute the sums Y^^Iq^^ \A"~^\ for all values of u in a linear (in 
n and d) number of operations. Thus, we only need a quadratic number 
of operations to find the z-th and the j-ih elements of A^j. 

Once x'^ and x^ are found, we can compute MaCy using Equation [14} 
We aim to find the coefficient of x^ in the polynomial (x^/x^)/,. Note 
that min{/ : x'^' divides x*} = min{/ : dj < ccj}. Therefore, it is easy to 
compute x'*/x^'. Now, the coefficient xl^ in (x*/x^')/,- equals the coefficient 
of x^/ (x*/x^') in //. This proves that we can find the coefficient we need 
in polynomial time. 

Since we have a polynomial-time algorithm that given n, d, i and / 
computes MaCy, there exists a polynomial-size circuit representing Mac. 

□ 

4.2 Determinant of a matrix given by a circuit 

In this subsection, we show that testing for zero the determinant of a 
matrix given by a circuit is PSPACE-complete. This result holds in any 
characteristic. The proof is based on the PSPACE-completeness of testing 
the existence of a unique s-f-path in a graph given by a circuit. 

The PSPACE-completeness of the s-f-connectivity (without the unicity 
condition) in a graph given by a circuit was proved via a reduction from the 
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problem QBF fyf[ . The idea was to turn a quantified formula into a graph. 
Since the proof of PSPACE-completeness of QBF actually uses a formula 
to express the s-f-connectivity in a graph, one may wonder whether their 
proof can be simplified by removing the use of a quantified formula. A 
by-product of our proof is a positive answer to this question. 

Lemma 2. Let C be a circuit representing a directed graph G with two distin- 
guished vertices s and t, with the promise that there is at most one path from s to 
t in G. It is PSPACE-complete to decide if such a path exists. 

Proof. The problem can be decided in polynomial space since the existence 
of a s-f-path in a (classical) digraph can be decided in nondeterministic 
logarithmic space. The classical algorithm is used as if the graph were 
given as input, and each time the algorithm needs to know if there is an arc 
from a vertex u to another vertex v, it evaluates the circuit on these inputs. 
This nondeterministic algorithm uses logarithmic space in the number of 
vertices, that is, polynomial space in the size of the input circuit. 

Now let L G PSPACE, decided by a deterministic Turing Machine M 
in polynomial space. Consider, on input x, the graph of configurations 
G^ of the machine M. This graph is described by a circuit that given as 
inputs two configurations c and c' of M, tells if c' can be reached from c 
in one step of computation. It is known that such a circuit of polynomial 
size exists, and that on input x, the circuit can be built is polynomial time. 
Now, there is a path from the start configuration to the accepting one in 
G^ if and only if x G L. 

To conclude, it remains to remark that there is at most one path from 
the start configuration to the accepting one in a deterministic machine. □ 

Corollary 2. Let C be a circuit representing a (directed) forest, that is a directed 
graph F such that each vertex has out-degree at most 1, and two distinguished 
vertices s and t. It is PSPACE-complete to decide if there is a path from s to t. 

Proof. The configuration graph of a deterministic machine is a forest. This 
shows that the previous proof actually implies this stronger statement. □ 

Let G be a digraph and M its adjacency matrix. A cycle cover of G is a 
subset of the arcs which form a set of cycles such that each vertex belongs 
to exactly one cycle. Therefore, if the vertices are numbered from 1 to n, 
a cycle cover corresponds to a permutation of {1, . . .,n}. The signature 
of the cycle cover is then defined as the signature of the corresponding 
permutation. It follows that the determinant of M equals the sum of the 
signatures of all the cycle covers of G. 

Theorem 7. Let C be a circuit describing a (0, \)-matrix M whose determinant 
is promised to he either 0, 1 or —1. Then it is PSPACE-complete to decide if this 
determinant vanishes. 
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Proof. A PS PACE algorithm for this problem follows from the (uniform) 
NC algorithm for the determinant [13]. 

Suppose we are given a circuit C describing a forest F, with two dis- 
tinguished vertices s and t. Consider the graph G obtained by adding an 
arc from f to s and adding loops on all remaining vertices. Then the only 
cycles in G are the loops and a cycle through s and t which exists if and 
only if there was a path from s to f in F. Therefore, a cycle cover of G 
exists if and only if there is a path from s to Hn F. 

If we consider the adjacency matrix of G, its determinant equals the 
sum of the signatures of the cycle covers of G. Then the determinant is 
nonzero if an only if there is a s-f-path in F, in which case its value is 
±1. □ 

Corollary 3. The problem of deciding if the determinant of a square (0, l)-matrix 
vanishes modulo n is PSPACE-complete for all n >2. 

Proof. For any n > 2, 1 ^ mod n and —1^0 mod n. Therefore, the 
theorem directly implies this corollary. □ 



5 Final remarks 

In characteristic zero, the upper and lower bounds on Resultant are in a 
sense close to each other. Indeed, NP = ZiP C AM C n2P, that is, AM lies 
between the first and the second level of the polynomial hierarchy. Fur- 
thermore, "under plausible complexity conjectures, AM = NP" [2, P157]. 
In particular, AM = BP-NPby definition, that is AM is a randomized ver- 
sion of N P, and randomization is often believed not to add any power to 
computation models. Improving the NP lower bound may be challenging 
as the proof of Proposition [1] shows that this would imply the same lower 
bound for Hilbert's Nullstellensatz. 

In positive characteristic, the situation is quite different. Indeed, the 
best known upper bound for Hilbert's Nullstellensatz as well as for the re- 
sultant is PS PACE. As in characteristic zero, the known upper and lower 
bounds are therefore the same for both problems. But as the gap between 
the NP lower bound and the PS PACE upper bound is rather big, these 
problems might be of widely different complexity (more precisely, test- 
ing the resultant for zero could in principle be much easier than deciding 
whether a general polynomial system is satisfiable). Canny's algorithm 
for computing the resultant [8] involves the computation of the determi- 
nants of exponential-size matrices, known as Macaulay matrices, in poly- 
nomial space. Those matrices admit a succinct representation. We proved 
that computing the determinant of a general succinctly represented matrix 
is PSPACE-hard (and testing it for zero is PSPACE-complete). It follows 
that Canny's polynomial-space upper bound could be improved only by 
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exploiting the specific structure of the Macaulay matrices in an essential 
way, or by finding an altogether different (non Macaulay-based) approach 
to this problem. As pointed out in Section ITT[ in characteristic zero a 
different approach is indeed possible for testing whether the resultant van- 
ishes (rather than for computing it). This problem is wide open in positive 
characteristic. 

Acknowledgments. We thank Bernard Mourrain and Maurice Rojas for 
sharing their insights on the complexity of the resultant in characteristic 0. 

References 

[i] E. Allender and M. Ogihara. Relationships among PL, #L, and the 
determinant. RAIRO-Inf. Theor. AppL, 3o(i):i-22, 1996. 

[2] S. Arora and B. Barak. Computational Complexity: A Modern Approach. 
Cambridge University Press, 1st edition, 2009. 

[3] J. Balcazar. The complexity of searching implicit graphs. Artificial 
Intelligence, 86(i):i7i-i88, 1996. 

[4] J. Balcazar, A. Lozano, and J. Toran. The complexity of algorithmic 
problems on succinct instances. Computer Science: Research and Appli- 
cations (R. Baeza-Yates and U. Manber, Eds.), pages 351-377, 1992. 

[5] R Biirgisser. Completeness and Reduction in Algebraic Complexity Theory. 
Algorithms and Computation in Mathematics. Springer, 2000. 

[6] P. Biirgisser, R Cucker, and R. De Naurois. The complexity of semi- 
linear problems in succinct representation. Computational Complexity, 
i5(3):i97-235/ 2006. 

[7] L. Buse and C. D' Andrea. On the irreducibility of multivariate subre- 
sultants. CR Math., 338(4):287-290, 2004. 

[8] J. R Canny. The complexity of robot motion planning, volume 1987 of 
ACM Doctoral Dissertation Award. MIT Rress, 1988. 

[9] J. R Canny. Some algebraic and geometric computations in PS PACE. 
In Proc. STOC'88, pages 460-469, 1988. 

[10] J. R Canny, E. Kaltofen, and L. Yagati. Solving systems of nonlin- 
ear polynomial equations faster. In Proc. SIGSAM'Sg, pages 121-128, 
1989. 

[11] J. R Canny and J. H. Reif. New lower bound techniques for robot 
motion planning problems. In Proc. FOCS'Sy, pages 49-60, 1987. 



22 



[i2] E. Cattani and A. Dickenstein. Introduction to residues and resul- 
tants. In A. Dickenstein and I. Emiris, editors. Solving polynomial 
equations, pages 1-61. Springer, 2005. 

[13] L. Csanky. Fast parallel matrix inversion algorithms. SIAM J. Comput., 
5:618, 1976. 

[14] C. D' Andrea and A. Dickenstein. Explicit formulas for the multivari- 
ate resultant. /. Pure Appl. Algebra, i64(i-2):59-86, 2001. 

[15] A. Dixon. The eliminant of three quantics in two independent vari- 
ables. Proc. Land. Math. Soc, 6:468-478, 1908. 

[16] I. Emiris and B. Mourrain. Matrices in elimination theory. /. Symb. 
Comput., 28(i-2):3-43, 1999. 

[17] J. Feigenbaum, S. Kannan, M. Vardi, and M. Viswanathan. Complex- 
ity of problems on graphs represented as OBDDs (extended abstract). 
In Proc. STACS'gS, page 216, 1998. 

[18] N. Fitchas, A. Galligo, and J. Morgenstern. Precise sequential and par- 
allel complexity bounds for quantifier elimination over algebraically 
closed fields. /. Pure Appl. Algebra, 67(i):i-i4, 1990. 

[19] H. Galperin and A. Wigderson. Succinct representations of graphs. 
Inform. Control, 56(3): 183-198, 1983. 



[20] M. R. Garey and D. S. Johnson. 
to the Theory of NP-Completeness. 
Sciences. W. H. Freeman, 1979. 



Computers and Intractability: A Guide 
Series of Books in the Mathematical 



[21] B. Grenet. Complexite du resultant et des grands determinants. Msc 
thesis, Ecole Normale Superieure de Lyon, 2009. 

[22] B. Grenet, P. Koiran, and N. Portier. The Multivariate Resultant is 

NP-hard in Any Characteristic. In P. Hlineny and A. Kucera, edi- 
tors, 35t/j International Conference on Mathematical Foundations of Com- 
puter Science, number 6281, pages 477-488, Berlin, Heidelberg, 2010. 
Springer- Verlag. 

[23] J. Heintz and J. Morgenstern. On the intrinsic complexity of elimina- 
tion theory. /. Complexity, 9:471-498, 1993. 

[24] D. lerardi. Quantifier elimination in the theory of an algebraically- 
closed field. In Proc. STOC'89, pages 138-147, 1989. 

[25] G. Jeronimo and J. Sabia. Computing multihomogeneous resultants 
using straight-line programs. /. Symb. Comput., 42(i-2):2i8-235, 2007. 



23 



[26] E. Kaltofen and P. Koiran. On the complexity of factoring bivariate 
supersparse (lacunary) polynomials. In Proc. ISSAC'o^, pages 208- 
215. ACM, 2005. 

[27] E. Kaltofen and P. Koiran. Expressing a fraction of two determinants 
as a determinant. In Proc. ISSAC'08, pages 141-146, 2008. 

[28] D. Kapur and T. Saxena. Comparison of various multivariate resultant 
formulations. In Proc. ISSAC'g^, pages 187-194, 1995. 

[29] D. Kapur, T. Saxena, and L. Yang. Algebraic and geometric reasoning 
using Dixon resultants. In Proc. ISSAC'94, pages 99-107, 1994. 

[30] P. Koiran. Hilbert's Nullstellensatz is in the polynomial hierarchy. /. 
Complexity, i2(4):273-286, 1996. 

[31] P. Koiran. Circuits versus trees in algebraic complexity. In Proc. 
STACS'oo, pages 35-54, 2000. 

[32] P. Koiran. The complexity of local dimensions for constructible sets. 

/. Complexity, i6(i):3ii-323, 2000. 

[33] P. Koiran and S. Perifel. VPS PACE and a transfer theorem over the 
reals. Computational Complexity, i8(4):55i-575, 2009. Preliminary ver- 
sion in Proc. STACS'07. 

[34] T. Krick, L. Pardo, and M. Sombra. Sharp estimates for the arithmetic 
Nullstellensatz. Duke Math. J., io9(3):52i-598, 2001. 

[35] S. Lang. Algebra. Springer, 3rd edition, 2002. 

[36] D. Lazard. Resolution des systemes d'equations algebriques. Theor. 
Comput. Set, i5(i):77 - 110, 1981. 

[37] A. Lozano and J. Balcazar. The complexity of graph problems for 
succinctly represented graphs. In Proc. 15th Int. Workshop WG'8g, page 
277, 1989. 

[38] F. Macaulay. Some formulae in elimination. Proc. Lond. Math. Soc, 
i(i):3, 1902. 

[39] F. S. Macaulay. The Algebraic Theory of Modular Systems. Cambridge 
Mathematical Library. Cambridge University Press, 1994. 

[40] G. Malod. Succinct Algebraic Branching Programs Characterizing 
Non-uniform Complexity Classes. In O. Owe, M. Steffen, and J. A. 
Telle, editors, Proc. FCT'ii, volume 6914 of Lecture Notes in Computer 
Science, pages 205-216. Springer, 2011. 



24 



[41] C. Papadimitriou and M. Yannakakis. A note on succinct representa- 
tions of graphs. Information and Control, 7i(3):i8i-i85, 1986. 

[42] D. A. Plaisted. New NP-hard and NP-complete polynomial and inte- 
ger divisibility problems. Theor. Comput. ScL, 3i(i-2):i25-i38, 1984. 

[43] J. Renegar. On the worst-case arithmetic complexity of approximating 
zeros of systems of polynomials. SIAM }. Comput., 18:350, 1989. 

[44] A. Seidenberg. A new decision method for elementary algebra. Ann. 
Math., 6o(2):365-374, 1954. 

[45] I. R. Shafarevich. Basic Algebraic Geometry 1: Varieties in Projective 
Space. Springer, 1994. 

[46] V. Shoup. New algorithms for finding irreducible polynomials over 
finite fields. NLath. Comput., 54(i89):435-447, 1990. 

[47] B. Sturmfels. Sparse elimination theory. In Proc. Comput. Algebr. Geom. 
Commut. Algebra. D. Eisenbud and L. Robbiano, eds., 1991. 

[48] J. Toran. Succinct representations of counting problems. In Proc. 
6th Int. Conf. Appl. Algebra, Algebraic Algo. and Error-Correcting Codes, 
pages 415-426, 1988. 

[49] B. L. van der Waerden. Modern Algebra, (third ed.) F. Ungar Publishing 
Co., New York, 1950. 

[50] J. von zur Gathen, M. Karpinski, and I. Shparlinski. Counting curves 
and their projections. Comput. Complex., 6(i):64-99, 1996. 

[51] K. Wagner. The complexity of combinatorial problems with succinct 
input representation. Acta Informatica, 2^{^):^2^-j56, 1986. 



25 



